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We report a precise measurement of the weak mixing angle from the ratio of neu- 
tral current to charged current inclusive cross-sections in deep-inelastic neutrino-nucleon 
scattering. The data were gathered at the CCFR neutrino detector in the Fermilab 
quadrupole-triplet neutrino beam, with neutrino energies up to 600 GeV. Using the 
on-shell definition, sin 2 # w = 1 - fj^V, we obtain sin 2 # w = 0.2218 ± 0.0025(stat.) ± 
0.0036(exp. syst.) ± 0.0040(model). 

PACS numbers: 13.15.Jr, 12.15.Mm, 14.80.Er 

The standard model (SM) of elementary particle physics describes the unification of 
the electromagnetic and weak interactions in terms of a weak mixing angle, sin 2 ^- In 
the on-shell convent ion QXJ] , the mixing angle is defined in terms of the W and Z boson 
masses: 

sin 9 W = 1 - — - T . (1) 

M z 

This Letter presents a SM extraction of sin 2 6\v from the ratio of neutral current (NC) to 
charged- current (CC) total cross-sections in deep-inelastic neutrino-nucleon (z/N) scatter- 
ing, 

+ nucleon — > + hadrons (NC), (2) 
Vy, + nucleon — ► fi~ + hadrons (CC). (3) 

The SM predicts that all electroweak processes may be described at lowest order in 
perturbation theory by just three independent experimental parameters. These may be 
chosen to be the electromagnetic fine structure constant (a), the Fermi coupling constant 
(Gf) and the mass of the Z boson (Mz), all of which have been measured to better than 
1 part in 10 4 . Electroweak processes cannot yet be predicted to this level of accuracy 
because higher order perturbative corrections for each process bring in additional depen- 
dence on the masses of the undiscovered top quark (M top ) and, to a lesser extent, the 
Higgs boson (Mni ggs )- Within the SM, the experimental determination of sin 2 #w from uN 
scattering has very little dependence on M top or MHi ggs @; in contrast, the SM prediction 
of sin 2 ^ from a, Gp and M z depends strongly on M top . Requiring the sin 2 #w from i^N 
scattering to agree with the prediction using M z sets limits on M top which are comparable 
with the best determinations from Z and W decay experiments at collider s[||. From a 
more general perspective, the consistency of the M top determinations from different pro- 
cesses constrains possible physics processes beyond the SM. Neutrino-nucleon scattering 
is uniquely sensitive to some proposed models with an extended Higgs sector or with 
extra Z's@]. Comparing the SM prediction for M\y from i/N scattering with the direct 
measurements at hadron colliders is a further test of the SM which is almost independent 
of M to p and M ffiggs . 



The E770 event sample of 3.1 million raw event triggers was collected in 1987-8 using 

w 



the quadrupole-triplet neutrino beam-line at Fermilab. The data sample for the sin 2 ^ 
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analysis contained 475627 events after all cuts, with a mean neutrino energy of 161 GeV 
and a mean 4-momentum transfer squared, Q 2 = 36 GeV 2 /c 2 . This represents approx- 
imately four times the statistics and almost twice the mean energy and Q 2 of the most 
precise previous sin 2 ^ determinations from i/N scattering^], 0|. 

The CCFR detector Q, [S| consists of a neutrino target /calorimeter followed by a muon 
spectrometer. The muon spectrometer was not used directly in the sin 2 #w analysis. The 
target comprises 168 iron plates, each 3 m x 3 m x 5.1 cm, interspersed with 84 liquid 
scintillation counters (every 10 cm of iron) and 42 drift chambers, each with x and y 
planes. It is 17.7m long, weighs 695 metric tons and has a mean density of 4.2 g/cm 3 . 

Both CC and NC interactions initiate a cascade of hadrons in the target that is 
registered by the drift chambers and scintillation counters. The muon produced in CC 
interactions typically penetrates well beyond the end of the hadron shower, appearing as 
a track of drift chamber hits with deposits of characteristic minimum-ionizing energies in 
the scintillation counters. We define the event length, L, to be the number of scintillation 
counters spanned by the event, where the longitudinal event vertex is defined to be the 
more upstream of the first 2 consecutive counters with more than 4 times the mean energy 
deposit of minimum ionizing muons ("mip's")0, and the event end is the counter above 
the next downstream gap of 3 counters with energies below 0.25 mip's. The mean position 
of the hits in the drift chambers immediately downstream from the vertex determines the 
transverse vertex coordinates. A calorimetric energy, E ca i, is calculated by summing up 
energy deposits in the 20 counters immediately downstream from the vertex. We require 
the event vertex to be more than 5 counters from the upstream end of the target and 34 
counters from the downstream end and less than 76.2 cm from the detector center-line. 
Requiring E ca i > 30 GeV ensures complete efficiency of the energy deposition trigger. 

The presence of a penetrating muon in CC interactions permits an approximate par- 
tition of CC and NC events by event length: 

^ _ No. events with L < 30 counters ^ No. NC 
No. events with L > 30 counters No. CC' 

where L > 30 counters implies a penetration greater than about 3.1m of iron. This 
experimental quantity was translated into a SM value for sin 2 #w using a detailed Monte 
Carlo-based computer simulation (MC) of the experiment which modeled the integrated 
neutrino fluxes, the relevant physics processes and the response of the CCFR detector. 
Predicted values for R meas were obtained by generating samples of simulated events and 
passing them through the same analysis procedure as the E770 data. The experimental 
value of sin^w was defined to be the input value to the MC which returned the same 
Rmeas as the E770 data. The relationship between R meas and sin 2 #w predicted by the MC 
is found, in a linear approximation, to be sin 2 ^ = 0.2218 - 1.73 (R^ c as - 0.4508). Our 
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experimental determination, R meas = 147795/327832 = 0.4508, corresponds to 

sin 2 fl w = 0.2218 ± 0.0025(stat.) ± 0.0036(exp. syst.) ± 0.0040(theor.). (5) 

The experimental and theoretical uncertainties were obtained from the MC by varying 
the model parameters within errors, and are itemized in Table 1. 

While the analysis presented here is performed completely within the context of the 
standard model, our value of sin 2 6w can be used with the MC model to calculate a cor- 
rected neutral to charged current event ratio corresponding to the incident z//z7flux[[|, R = 
(NC U + NC 17 )/(CC U + CC 17 ) = 0.3117±0.0014(stat.)±0.0018(exp. syst.) ±0.0014(theor.). 
This value corresponds to a hadron energy cut of 30 GeV on both CC and NC events. The 
event ratio is fully corrected for experimental effects such as acceptance, smearing, and the 
v e background but no theoretical corrections are applied other than an isoscaler correction. 
For the quantity R, the theoretical uncertainty is almost entirely due to the longitudinal 
structure function, R long . In terms of R v & = a^}/^cc\ R ~ 0.895^" + 0.105R 17 . The 
variation of R with sin 2 9w is very similar to that of R meas and is given by dR/d sin 2 0w = 
-0.565. 



The integrated and fluxes for E770 were determined directly from low hadron 
energy CC event samples, normalized to the neutrino total cross sections [fDJ. The final 



event sample consisted of 86.4% v^, 11.3% and 2.3% v c or z/ e interactions. Errors 
in the flux tend to cancel in the ratio R m eas, but the z/ e flux modeling is more critical 
because essentially all z/ e events are short enough to appear in the numerator of R m eas- The 
integrated z/ e flux was modeled using a Monte Carlo simulation of the neutrino beam-line, 
with the spectra of secondaries from the proton target parameterized from experimental 
production cross- sections fTTH . Approximately 80% of the z/ e 's in the final data sample 



were produced from the Ke3 decay mode of charged kaons, whose modeling is directly 
related to the observed event spectrum. The next largest contribution to the z/ e flux 
was from neutral kaon decays (~16%), with smaller contributions from the decays of D 
mesons, pions, muons, A's and E~'s. 

The modeling of neutrino-induced events in the detector and resolution smearing ef- 
fects on the measured L, E ca i and vertex positions were modeled primarily using neutrino 
and test beam data events[[F], |8], [T2| . Systematic uncertainties associated with the hadronic 



energy measurement E ca i include possible small NC/CC shower differences, uncertainties 
in the muon energy deposit within the hadron shower, and uncertainties in the resolu- 
tion function, e/7r response, and absolute energy scales obtained from hadron/electron 
test beam measurements [^, The length uncertainties include those associated with 



the shower length parameterizations from test beam measurements p2| , the longitudi- 
nal vertex determination, which has been checked against the vertex of dimuon events, 
counter inefficiencies and noise, and effects from dimuon production. Small differences 
in the transverse vertex for NC and CC events due to the muon drift chamber hits were 
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determined from the analysis of CC events with these hits removed. 



The NC and CC differential cross-sections were modeled using a QCD-enhanced quark- 
parton description of the nucleon. The quark distributions were obtained by using a 
modified Buras-Gaemers parameterization Jl3| of CC nucleon structure functions measured 



in the same CCFR experiment ||14|| . The strange quark component, parameterized by the 
momentum fraction relative to the non-strange sea, k = 2s/(u + d), was determined from 
an analysis of CCFR dimuon events [T^] which arise from the muonic decays of charm 
quarks produced in CC scattering off down and strange quarks and anti-quarks. The 
threshold suppression for this process, due to the mass of the charm quark, was modeled 
using a leading order slow-rescaling formalism, with a fitted effective charm quark mass of 
M c = 1.31 ± 0.24GeV/c 0]. The level of the charm 

sea was assumed to be 10% of the 



strange sea, consistent with a wrong-sign muon analysis from a previous CCFR neutrino 
experiment using the same detector [[U|. Our parameterization of the Ri ong is based on 



QCD predictions and data from charged lepton scattering experiments ||17|| . A correction 
for the difference between u and d valence quark distributions in nucleons, obtained from 
muon scattering dataflSfl , was applied to account for the 5.67% excess of neutrons over 



protons in the target. Radiative corrections to the scattering cross-sections were applied 



using computer code supplied by Bardinjnj, assuming values of M top = 150 GeV/c and 
M ffi ggs = 100 GeV/c 2 . 

Figure 1 shows the length distribution of the E770 final data sample and a MC simu- 
lated event sample. Events reaching the muon spectrometer, comprising 79% of the CC 
interactions, have been left out for clarity but are included in the normalization of the 
MC event sample to the data. The remaining CC events have a muon which either has 
a low energy and ranges out in the neutrino target or has a large opening angle with 
respect to the incident neutrino and exits through the side of the target. The produc- 
tion energy and angular distributions of these muons are very well constrained by the 
CCFR structure function measurements, and their propagation through the target has 
been precisely parameterized using large samples of muons from test beam and neutrino 
data. The events with length less than or equal to the 30 counter partition of equation 
4 are predominantly true (or V^) NC events, with 22.9% and 7.3% backgrounds from 
short and CC events and u e events, respectively. Since the NC and u e event lengths 
fall well short of the 30 counter partition, the sensitivity to the modeling of the hadron 
shower length is minimal. The good agreement between data and MC for event lengths 
greater than 25 counters reinforces confidence in the estimate of the CC component of 
the event sample in the 'NC length region of 30 counters or less. In addition, the data 
and MC distributions in E ca i and vertex radial position also agree well. 

The most precise previous determinations of sin 2 #w in i/N are from the CDHS[|| and 
CHARM collaborations. After adjusting to our theoretical assumptions for the charm 
quark mass, M c = 1.31 ± 0.24 GeV/c 2 |15| , and top quark mass, M top = 150 GeV/c 2 , 
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these experiments yield measurements of sin 2 6\y = 0.2225 ±0.0066 (CDHS) and sin 2 6\y = 
0.2319 ± 0.0065 (CHARM), in agreement with our result, sin 2 # w = 0.2218 ± 0.0059. 

Combining our value of sin 2 #w with the precise measurement of the Z boson mass, 
M z = 91.187±0.007GeV/c 2 @, gives M w = 80.44±0.31 GeV/c 2 , and corresponds in the 
SM to M top = 1901 

48^14 GeV/c 2 pi] , where the central value and first set of uncertainties 
assume Mmggs = 300 GeV/c 2 and the second set of uncertainties come from varying Mmggs 
between 60 and 1000 GeV/c 2 . Our results are consistent with the values, Mw = 80.25±0.10 
GeV/c 2 and M top = imt{g±H GeV/c 2 , from a SM fit to a large number of experimental 
results from Z decays at the LEP electron collider |J, and our Mw value also agrees with 



the world average measured value from hadron colliders, Mw = 80.14 ± 0.27 GeV/c 2 [22 



In summary, our value is the most precise determination of the on-shell sin 2 #w i n a 
single experiment and is consistent with previous determinations in uN scattering and 
other processes. 
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Foundation and the Department of Energy. SRM acknowledges the support of the Alfred 
P. Sloan Foundation. 
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TABLE 1. Uncertainties in the measurement of sin^w- 



data statistics 


0.0024 


Monte Carlo statistics 


0.0006 


TOTAL STATISTICS 


0025 


muon neutrino flux 


U.UUUO 


[y e ± 4.2% ) electron neutrino flux 


U.UUzo 


transverse vertex 


u.uuuy 


(±25%) cosmic ray subtraction 


U.UUUo 


Energy Measurement 




NC/CC shower difference 


n nnn7 

U.UUU 1 


muon energy loss in shower region 


U.UUUO 


IZC1U/0J llctLli Ull C11C1 g V 1 coUlll L1U11 


U.UUUO 


(±1%) absolute energy scale 


U.U018 


Event Length 




hadron shower length 


0.0010 


(±2.5 cm) vertex determination 


0.0010 


counter efficiency and noise 


0.0009 


dimnori nrodnction 


0003 


TOTAL EXP SYST 


0036 


(M c = i.ol ± U./4 LieV/c j cnarm proa. 


U.UUoU 


< R, 7C 1^% lr.no- QTT 

I ^o«9+ J-o/o J long, 


n nm q 
u.uuiy 


(C/S = 0.10 ±0.15) charm sea 


0.0015 


rad. corrections 


0.0007 


higher twist 


0.0005 


non-isoscalar target 


0.0004 


(k = 0.37 =F 0.05) strange sea 


0.0003 


structure functions 


0.0003 


TOTAL PHYSICS MODEL 


0.0040 
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Figure Captions 



FIGURE 1. Data and Monte Carlo (MC) event length distributions. The data are 
represented by dots and the MC prediction by the solid line. Also shown are the MC 
contributions from NC events ("NC"), CC events ("CC") and combined NC and 
CC interactions from v c or V c (V e "). The inset shows the data, total MC and the NC 
contribution to the MC for the region L>25 counters. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ex/9405008vl 



